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Searching for a marked target in a list of N items requires Sl(N) oracle queries when using a 
classical computer, while a quantum computer can accomplish the same task in 0(y~N) queries 
using Grover's quantum search algorithm. Although this quadratic speed-up has been indisputable, 
the Copenhagen interpretation of quantum mechanics does not make understanding how this en- 
hancement over classical occurs intuitive. At the same time, cosmologists have proposed the theory 
of decoherent histories as an alternative to the Copenhagen school in order to explain how classical 
physics arises from the more fundamental quantum laws. This article applies theory of decoherent 
histories to Grover's algorithm. It finds a quasi-classical domain within the search evolution and 
compares it to classical searching. Moreover, under the theory of decoherent histories, the quantum 
process can be described in a language very much similar to a stochastic process. Thus, quan- 
tum searching, which has been of primal importance in quantum computing, also gives an exactly 
£SJ ■ solvable toy model for the theory of decoherent histories. 
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PACS numbers: 03.65.Sq, 03.67.Ac 



d . I. INTRODUCTION. 

Quantum computer science has come to know numerous examples in which a quantum computer can outperform 
its classical counterpart [l[. Among the celebrated algorithms is that a quantum computer can search for a marked 
element from a list of N items using the oracle only 0(\/~N) number of times rather than the classical il(N) 0. 
The quantum algorithm is entirely dissimilar to the simple process of elimination of classical physics. We have an 
■ initial state of equal superposition of all candidates that rotates into the target after a time proportional to 0(sfN). 
CTN \ This picture is rather unintuitive, and a deeper understanding of how the quadratic speed-up happens is needed. 
t-H . Entanglement has often been considered the key feature responsible for quantum-over-classical enhancements [H, 0] , 
\q j but in the case of searching, it has been shown special purpose devices can enjoy the quadratic speed-up without 
entanglement [f|. However, it has also been argued interference or quantum correlation may be the key resource 0,0|- 
Thus, to identify the essential non-classical features responsible for the quantum enhancement remains a challenge. 

On the other hand the division between classical and quantum can be a flaw in interpretation. Copenhagen quantum 
mechanics posits two kinds of time evolution: unitary and collapse of the wave function when the system interacts 
with the measuring apparatus. It is as if the world were composed of quantum and classical devices. However, we 
know even our classical computers are fundamentally quantum. Such a distinction has been found unsatisfactory in 
both logic and aesthetics by people interested in foundational issues of quantum mechanics as well as cosmologists, 
who have invented the theory of decoherent histories in hopes of explaining how classical physics arises from quantum 
0, Q. In this paradigm the measuring apparatus is included with the system of interest so as to form a closed 
quantum system. The goal is to assign probabilities to events and to find correlations in time of an observable. If a 
certain evolution path of that observable is overwhelmingly more likely than any others, then it is said to follow a 
deterministic classical law. For example, among the numerous correlations of positions of the moon at different times, 
the one that is most likely should be that which is consistent with Newton's law. Of course, such a program is so 
ambitious that it has not yet been done, if ever, and the literature has but a few examples illustrating how the theory 
works dillD]. 

Since the theory of decoherent histories has the potential to derive classicality from quantum, it has been applied 
to quantum information processes in hopes of discovering how quantum computing differs from classical computing. 
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In [12j Poulin investigated extending quantum algorithms, which typically consist of the initial state preparation, the 
dynamics, and the final measurement, to multi-time families of decoherent histories with intermediate measurements. 
It demonstrated through a toy example how the extension is done and provides a number of theorems concerning the 
general structures of such multi-time extensions. In the current paper we apply the theory of decoherent histories 
to Grover's algorithm as a particular example illustrating how a classical process can arise within a closed quantum 
system. We then compare it to the classical search. Moreover, we are able to compute the main figure of merit in 
closed form, thus making the Grover's algorithm an exactly solvable toy model for the theory. Studying the quantum 
search in this paradigm allows one to describe the process as a motion picture in the language similar to a stochastic 
process. The beginning is very much a classical deterministic evolution in which certain meaningful events can be 
assigned probabilities. Then, interference between alternative histories starts to become significant that we are in the 
"quantum" regime, after which the target is found. In the next section we will review the basics of the theory so as 
to set the stage before presenting our analysis with the main findings. Finally, we will conclude with some possible 
future directions of research. 



II. DECOHERENT HISTORIES 



One of the primary goals of the theory of decoherent histories is to recover and demonstrate how classical physics 
arises from the more fundamental quantum laws. Thus, let us study a closed quantum system, and if there is a 
measuring apparatus, we would include it with the system of interest and study the whole composite. Unlike the 
Copenhagen interpretation, the theory of decoherent histories does not assume the measuring apparatus to be classical 
whose dynamics collapses the wave function (see [13l - ll5| for a complete review). Suppose the system evolves according 
to the unitary evolution U(t) and starts in the state \s). Consider events at different times formed by a resolution of 
the identity. For example, let Pi with i £ {1, . . . , M}, be orthogonal projectors such that 

M 
i=l 

P t P 3 = (2) 

The projectors serve as an algebra of events where union and intersection are achieved by addition and multiplication, 
respectively. At a different time we may consider a different resolution of the identity and thereby speak about a 
different observable. The time intervals between different resolutions of the identity can be chosen arbitrarily. A 
family of histories consists of choosing such sets of orthogonal projectors for different times up to a certain total time, 
and a history is a tensor product of projectors at different times. For the subsequent analysis let us consider the same 
resolution for all times and the same amount of time, t, between successive events for a total number of K times. 
Denoting i — (ii, %2, ■ ■ ■ , %k) € {1, ■ • ■ , M} K , we define an elementary history to be 

E- := P lK ® P lK _ 1 ® • • • ® P it . (3) 

The history above has the meaning that at time t, Pi t happens, at time 2t, Pi 2 happens, and so on and so forth such 
that at time Kt 1 Pi K happens. 

Having defined a family of histories, we would like to be able to assign probabilities to them. The general state of 
affairs is that it is not always possible to do so, and even when it is, we have a stochastic process which may or may 
not be classical. In the Young double-slit experiment, for example, we have a family of 2-time histories with the first 
being concerned with the event which slit the particle goes through and the second the particle's position. It is well 
known that if we could assign definite probabilities to the first event, probability sum rules would go horribly awry. 
On the other hand, a case where we can assign probabilities is Grover's quantum search algorithm. Consider a family 
of 2-time histories based on the events whether we find the target at exactly the times when the state vector is the 
target. Trivially, we have one history of unity probability, namely that in which we find the target both times, and 
the rest zero. Although it makes perfect sense to speak of probabilities for different histories in this family, we would 
not consider such a process classical. After all, it takes only 0{y~N) queries. Thus, being able to assign probabilities 
to a family of histories is necessary but not sufficient for an evolution to be regarded classical. 

The reason probabilities can not be assigned is that the usual probability calculus may fail. Recall that we expect 
the probability of a union of exclusive events to be the sum of the probabilities of the individual events. Given a 
family of histories we can have unions of exclusive histories by forming a coarser- grained family consisting of sums of 
elementary projectors belonging to certain time slices. For example, at time kt we may form the following resolution 
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with je {!,..., M-l}, 

P> = { P i Sj<M-l, 
] \Pm-i+Pm Hj = M-l, K ' 

The histories in this new family differ from those of the old one in their fc-th component. An example history might 
thus be 

F = P iK ® • • • ® (P M -i + Pit) <8 ■ ■ ■ ® Ph , (5) 
which is the union of the two exclusive histories of the old family 

\e ? =P iK ®---®P M ®'--®Pi L . 

Thus, if we have a rule to assign probabilities to the histories E*, then we expect 

Probability^) = Probability^-) + Probability (Ep) (7) 

consistent with the usual sum rule for disjoint events. We can form numerous other coarser-grained families by 
summing the elementary projectors at additional other time slices, and in all of them, the probability of a coarser- 
grained history should be the sum of probabilities of the finer-grained histories. 

The necessary and sufficient condition for probabilities to be assigned successfully has been identified 14j. Corre- 
sponding to each history i, we define the following operator 

C I :=P lK U(t)P lK _ 1 U(t)---P ll U(t), (8) 

and the decoherence functional D be defined as 

D{li>) {s\ClC-\s) (9) 

for any pair of histories i and V . It is possible to assign probabilities to histories if and only if for all i ^ V 

Re(D(i,^)) = 0. (10) 

Moreover, the probability of history i would be D(i,i), which has also been demonstrated to be consistent with Born's 
rule. By Equations © and ©, the probability of the coarser-grained history F would be 

Pr(F) =Pr(%)+Pr(£7) + 2Re(Z)(T,i 7 ')). (11) 

Thus, Re(D(i, V)) — allows us to recover the usual probability calculus. 

The importance of the exact decoherence condition, Equation (|10p . is that the family giving rise to it is a stochastic 
process. If we are to show classical behaviour emerging from a closed quantum system, we necessarily need to find 
a family of histories satisfying Equation (|10[) . It would also be useful to consider the less stringent approximate 
decoherence condition that for all i =/= i' 

Re(D(i,i?))xtO, (12) 

in which case probability sum rules would be violated slightly. Precisely how small is small enough for probabilities to 
be assigned meaningfully depends on the situation [16j . However, it has been shown [13l | if there exists an < e < 1 
such that for all i ^ V we have 

Re(D{i,i ! )) 



y 'D(i,i)D(i> ',?) 



< e, (13) 



then probability sum rules for most coarse grainings are satisfied to order e. To understand precisely what 'most' 
means, one would need to refer to [l3j], but the sense in which probability calculus is obeyed is that, for our simple 
example of coarse graining above, 

Pr(F) < (1 + e)(Pr(%) + Pr(£? i7 )). (14) 

Approximate decoherence has been shown to arise for local densities connected to globally conserved quantities [l8j . 
Indeed, we will see in the next section that the Grover's search algorithm admits a globally conserved quantity, and 
a local density gives rise to a quasi-classical regime. Moreover, the proposed measure of classicality in Equation (|13[) 
will be shown to be too strict and can fail unnecessarily in some cases. 
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III. ANALYZING GROVER'S QUANTUM SEARCH ALGORITHM 

Let us reformulate the continuous time version of Grover's quantum search algorithm (l7j so as to give it a physical 
context that is amendable to analyzing it using the theory of decoherent histories. We will study the time correlations 
of a local density of a globally conserved observable and see that its evolution is quasi-classical. Consider the N 
dimensional Hilbert space of N independent spins spanned by {|m)}^ =1 where |m) denotes spin m being up while 
the rest of them down. The marked element w which is to be found is distinguished in the following manner in the 
Hamiltonian 

N 

H w = -J2(-l) 5 ™a m -(N-2) ( 15 ) 
m=l 

where a m , in the ordered basis is the tensor product of ^ on the m-th spin and identity on the rest. 

This Hamiltonian H w has the desired spectrum of all the ground states having energy expect for the marked state 
having energy 2. To this Hamiltonian we add the driving Hamiltonian 2 |s)(s|, where 

1 N 

\s) = — J2\m) (16) 

v in— 1 

and consider henceforth the total Hamiltonian 

H = H w + 2\s)(s\. (17) 

The system starts in the state \s) and evolves unitarily according to the Hamiltonian H. As was shown in [17j . this 
system will find the unknown target wina time 0(y/~N). 

To find a classical stochastic process within this system, we need to see for which families of histories the decoherence 
functional D becomes vanishingly small. It has been shown that a quasi-classical domain comes from local densities 
connected to a globally conserved quantity [l8[ . In our case the globally conserved quantity is the total magnetization 
Sm=i °Vni an d we are thus led to consider local densities. Although there are numerous choices for a local density, 
we choose to consider J2 m ^w a m since its eigenspace projectors 

[p =HH (18) 
1 p = i-p y ' 

represent the events of success finding the target and failure, respectively. Using this resolution of the identity, we 
will build a family of histories with time t between successive events for a total of K times. We will take t to be 
equivalent to 1 application of the oracle in discrete time and K of the order 0(-\fN). For brevity let us denote success 
by 1 and failure by 0. Then, a history is represented by a € {0, 1} K ■ When successive values are alike, we can collect 
them into a streak. For example, if a history consists of 0's for the first si events, l's for the next S2 events, and so 
forth up to l's for the last s n events, then we would write it as a = I s ™ . . . l S2 Sl , where n is the number of streaks 
in the history with 1 < n < K and Sj the length of streak j. Hence, 



K. (19) 



Exact results 



With the operators C of Equation © we have the unitary evolution 

U{tK)\ S ) = {{P + P)U{t)) K \s) (20) 
" C a \s) . (21) 



o£{0,l} K 



Instead of considering the unitary evolution U(tK) on the left hand side of Equation (|2U|) as is traditionally done, we 
can consider the same process as a sum over all possible histories a as shown on the right hand side of Equation (|21l) . 
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At each time step, after one application of the oracle, we ask the question: "Is the target found?" Each history a is a 
string of possible answers "yes" or "no" to the question and, thus, is very much like a path in a random walk except 
that we may not be able to assign traditional probabilities to the paths. As discussed in the previous section, we need 
to see if the decoherence functional D vanishes for any two distinct histories in this family. To this end, let us define 
the following quantities in terms of N and w. 
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1 



E 



\m) 



m={l,...,N}\{w} 

f := cosxt — 12; sin art = l/le 1 ^, 
One can then show that, with much algebra skipped, we have 



1 



N 



N-l . 
b := — — — sinri. 



(22) 
(23) 



C a \s) = \f\ K - n b n -U 



exp(i<j>(^ ^2 a k- E Qfc )) (/ V N n 1 + bx ) ^ ' if a = I s ™, • • ■ ,0 si ,n = even 

fc-cven A;— odd 

exp(i^( a k - E a k))(b^^^+\f\x)\w), ifa=l s »,...,l*i,n = odd 

k— odd k— even 



explKpl 2^ a k 

k— odd k— even 

ex p( i( ^( E afe_ E afc ))(i/i 

fc-even k— odd 



N- 1 

N 



if a = s " , . . . , l Sl , n — even 



+ e I0 &x) if a = O s ",...,0 Sl ,n = odd 



We thus have derived an exact expression for the different histories of the analog version of Grover's quantum search 
algorithm. Next we will determine how the local spin density ^2 m ^ w o m exhibits a quasi-classical evolution. 



B. Quasi-classical domain 

From here on let us work in the regime N 1. It is helpful to recognize the magnitudes of the different relevant 
quantities. Since x — 1/yN and t is the duration equivalent to one application of the oracle in discrete time, we have 
xt e O(x) = 0(1/ VN) = o(l). Thus, 



|6| w xt 
|0| » 1/JV 



(25) 



Moreover, the total time for which we study the process is proportional to K, so long as K < V N, the phases are of 
order <pK ps 1/\N and therefore negligible. As a result, by Equation (|24|). we get 



l/l 



K— n lti—1 



' (/* + bx} \w) 
b + \f\x)\w) 
b + fx) |0, 
\.f\ + bx] |0 , 



if a= I s ", 
if a= I s ", 
if a = s - , 
if a = s " , 



. , S1 , n = even 
. , l Sl , n — odd 
. ,l 3l ,n = even 
. ,0 Sl ,n = odd 



(26) 



In the beginning when K = 1, we have only 2 branches with the dominant being of amplitude approximately 
unity and the other 1 being approximately zero. Even when K increases to a few, Q K is still the most dominant as 
other branches come into existence of amplitude at most 0(b). In other words, for a short time in the beginning 
the local density J2 m ^w ° m mos t surely has value — N + 3 with small quantum fluctuations. This regime is highly 
deterministic and should be considered quasi-classical. However, when we compute the measure of classicality in 
Equation (fT3|) between the small branches, we would get e ps 1 due to its being a ratio of two small yet comparable 
quantities. Although probability sum rules would be violated to a great extent, we nevertheless consider this regime 
classical since the cases for which violations occur are extremely unlikely. Thus, this measure of Equation (|13p does 
not quite reflect our intuition that the system should be considered quasi-classical at this stage. 
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C. Time of Grover's Search Algorithm 

As K increases the number of different histories a grows exponentially and the amplitude of each branch decays 
exponentially so as to keep normality. Eventually 0^ gets killed by time and becomes comparable to other branches. 
Let us now consider the probability of success at a given time. As shown by Equation (|21[) it comes from the sum of 
branches that end in 1, which can be broken up as 



A" 



E cr- = E E 



C a r, 



(27) 



histories of n streaks 
(K-l\ 



For a given K, the number of histories with n streaks is Thus, with the results of Equation (|26p . the sum of 

branches that end in success is (with approximations taken into account in the limit of iV > 1) 



E c « I s ) 

ae{l}x{0,l} A ' _1 



where 



K 

E 



K - 1 
n — 1 



|/|*-*6"- 1 5(n) \w) 



(28) 



S(n) = 



f* , n = even 

b + \ f\z, n = odd, 



(29) 



while the sum that ends in failure is 



where 



K 



E cr«W = E 

a£{0}x{0,l} K - 1 



n=l 



K- 1 
n-1 



\f\ 



K—n j n-l 



F(n) \0 



(30) 



F(n) 



b + fx, n = even 
\f \ + bx, n = odd. 



(31) 



When K G O(l) the value n that maximizes ( K ~^\\f\ K ~ n b n ~ l is near 1, thus making the branch that ends in failure 
dominant. As K increases, the critical n moves away from 1 and when it hits 2, the sum of branches that end in 
success can be approximated by 



E c « 



K — 1 
1 



\f\ K - 2 b(f*), 



(32) 



due to sub-dominant terms being dropped out, while that which ends in failure is 



(£iE c ° 



K — 1 
1 



\f\ K - J b(b + fx). 



(33) 



Hence, at this stage the probability of success is nearly 1 because b + fx G 0(1/ VN). See Figure[T]for an illustration 
of the relative sizes of the terms in the summations of Equations (l28l) and (l30l) at the time when the target gets found. 

We are thus interested in finding the critical n that maximizes (^Ti) l/| X_ra |^| n_1 m the regime n G 0(1) and 
K ^> 1 as a function of K . Using Stirling's approximation, we find the critical n to satisfy 



logl^-^Q 
S |/|(n-l) 

Since we are dropping sub-dominant terms, the time at which the target is found is approximately 

\b\ 

which recovers the well-known time that Grover's search algorithm requires to succeed. 



(34) 



(35) 
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FIG. 1: Plotted are the terms in the sums of Equations (|28[) and (|30f) for a database size of N — 100 and K = 9, roughly 
corresponding to the time when the target is found. The y-axis is the amplitude of a history and the x-axis the number of 
histories of a certain number of streaks. The sizes of the bubbles correspond to the numbers of streaks with the smallest being 
1 and the largest 9. The green bubbles are histories that end in success, red in failure, and black as reference points whose 
product of x and y coordinates is 1. The product of a bubble's x and y coordinates is the magnitude of its contribution to 
the sum. Histories that oscillate between failure and success are suppressed either in amplitude, the quantity, or both. In fact 
only histories with 3 or fewer streaks contribute to the sums substantially. The most dominant history is a streak of failures 
followed by success, which explains why at this time the target is found. 



IV. CONCLUSION 

Wc find studying the Grover's algorithm in the framework of decoherent histories allows us to describe the search in a 
language very much like a stochastic process. Moreover, it provides a nice and exactly solvable toy model of the theory 
in which one can see how a quasi-classical regime arises out of the more fundamental quantum laws. The dynamics 
of Grover's quantum search algorithm admits a local density whose evolution consists of one dominant history with a 
few much smaller alternatives. During this quasi-classical regime the probability assigned to the dominant history is 
nearly 1 and is the same as the classical search to the degree of accuracy computed. As time goes on, histories of all 
different streaks of O's and l's begin to populate-in fact, this number grows exponentially with time-until histories 
of the type a streak of O's followed by a streak of l's become dominant. At which point the target is found. Thus, 
the theory of decoherent histories provides for a motion picture of how Grover's search algorithm finds its target in a 
language similar to a stochastic process. Perhaps it is more intuitive than thinking about the process as a rotation of 
an abstract state vector and hopefully contributes to further work on pinpointing exactly the key feature of quantum 
theory that is responsible for the quadratic improvement. 
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